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Now that the properties of low temperature Bose gases at low density, p, can be examined experi- 
mentally it is appropriate to revisit some of the formulas deduced by many authors 4-5 decades ago. 
One of these is that the leading term in the energy /particle is pa/m, where a is the scattering 
length. Owing to the delicate and peculiar nature of bosonic correlations, four decades of research 
have failed to establish this plausible formula rigorously. The only known lower bound for the energy 
was found by Dyson in 1957, but it was 14 times too small. The correct bound is proved here. 
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With the renewed experimental interest in low density, 
low temperature Bose gases, some of the formulas posited 
four and five decades ago have been dusted off and re- 
examined. One of these is the leading term in the ground 
state energy. In the limit of small particle density, p, 



eoip) « fi An pa, 



(1) 



where eo(p) is the ground state energy (g.s.e.) per parti- 
cle in the thermodynamic limit, a is the scattering length 
(assumed positive) of the two-body potential v for bosons 
of mass TO, and p, = h?' /2m. 

Is (^ correct? In particular, is it true for the hard- 
sphere gas? While there have been many attempts at a 
rigorous proof of (|^) in the past forty years, none has been 
found so far. Our aim here is to supply that proof for 
finite range, positive potentials. As remarked below, (|^) 
cannot hold unrestrictedly; more than a > is needed. 

An upper bound for eo(p) agreeing with (|l|) is not 
easy to derive, but it was achieved for hard spheres by 
a variational calculation which can be extended to 
include general, positive potentials of finite range. What 
remained unknown was a good lower bound. The only 
one available is Dyson's 0], and that is about fourteen 
times smaller than (Q). In this paper we shall provide 
a lower bound of the desired form, and thus prove (0). 
We can also give explicit error bounds for small enough 
values of the dimensionless parameter Y = 47rpa'^/3: 



eo{p) > p Anpail - CY^^") 



for some fixed C (which is not evaluated explicitly be- 
cause C and the exponent 1/17 are only of academic in- 
terest). The bound (|^) holds for all non-negative, finite 
range, spherical, two-body potentials. A typical experi- 
mental value M is F « 10^^. Dyson's upper bound is 
p 47rpa(l + 21^/3)(1 - yi/^)-^. 

We conjecture that (^ requires only a positive scatter- 
ing length and the absence of any many-body, negative 
energy bound state. If there are such bound states then 
(^ is certainly wrong, but this obvious caveat does not 
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seem to have been clearly emphasized before. There is a 
'nice' potential with positive scattering length, no 2-body 
bound state, but with a 3-body bound state [1). 

Our method also obviously applies to the positive tem- 
perature free energy (because Neumann boundary con- 
ditions give an upper bound to the solution to the heat 
(or Bloch) equation). 

We also give some explicit bounds for finite systems, 
which might be useful for experiments with traps, but we 
concentrate here on the thermodynamic limit for simplic- 
ity. For traps with slowly varying confining potentials, 
Vcxt, our method will prove that the leading term in the 
energy is given by the well known local density approxi- 
mation [|, which minimizes the gaseous energy (|l|) plus 
the confining energy, with respect to p(x), namely. 

Sip) = J (V;,t(x)p(x) +Ai4^ap(x)2)d3x 

is minimized subject to J p — N — number of particles. 

The fact that Dyson's lower bound was not improved 
for four decades, despite many attempts, attests to the 
fact that bosons are subtle quantum mechanical ob- 
jects which can have peculiar correlations unknown to 
fermions. For example, there is the non-thermodynamic 
iV^/^ law for the charged Bose gas that was discovered by 
Dyson confirmed only 20 years later 1^, and which 
defies any simple physical interpretation. 

The first understanding of (|l|) goes back to Bogoliubov 
who also introduced the notion of 'pairing' in Helium 
(which resurfaced in the BCS theory for fermions). Later, 
there were several derivations of (^ (and higher order) 
|9j. The method of the pseudopotential, which is an 
old idea of Fermi's, was closest to the Bogoliubov analy- 
sis. The 'exact' pseudopotential was constructed in pO[ , 
but it did not help to make this appealing idea more rig- 
orous. Most of the derivations were in momentum space, 
the exception being joj, which works directly in physical 
space and which can handle both long and short range 
potentials. See [|ll[ for a review. All these methods rely 
on special assumptions about the ground state (e.g., se- 
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lecting special terms in a perturbation expansion, which 
hkely diverges) and it is important to derive a fundamen- 
tal result like (P without extra assumptions. 

In all of this earlier work one key fact was not under- 
stood, or at least not clearly stated in connection with 
the derivation of It is that there are two different 
regimes, even at low density, with very different physics, 
even though the simple formula (|l|) seems to depend only 
on the scattering length. Recall that the (two-body) scat- 
tering length is defined, for a spherically symmetric po- 
tential, D, by 

-fiu'i{r) + ^vir)uo{r)=0, (3) 

with uq(0) — 0, uo{r) > (which is equivalent to the 
absence of negative energy bound states, and which is 
true for nonnegative v). As r — > oo, u{r) « r — a. (Note 
the v/2 and not v in (^ because of the reduced mass.) 
Thus, a depends on m in a nontrivial way, and there are 
two extremes: 

1. Potential Energy Dominated Region: The hard 
sphere (v(r) = oo for r < a), is the extreme case here; 
the scattering length is independent of m, and the energy 
is mostly (entirely) kinetic. We see this from (P because 
—m deo/dm is the kinetic energy (Hellmann-Feynman 
theorem). In this regime the potential is so dominant 
that it forces the energy to be mostly kinetic. The g.s. 
wave function is highly correlated. 

2. Kinetic Energy Dominated Region: The typical case 
is a very 'soft' potential. Then a « (m/h^) v{r)r'^dr, 
which implies, from (|^), that cq hardly depends on m. 
Thus, the energy is almost all potential. The g.s. wave 
function is essentially the noninteracting one in this limit. 

In other words, 'scattering length' is not a property of 
V alone, and the low density gas, viewed from the per- 
spective of the bosons, looks quite different in the two 
regimes. Nevertheless, as (|l|) says, the energy cannot dis- 
tinguish the two cases. Whether Bose-Einstein conden- 
sation itself can notice the difference remains to be seen. 
Condensation will not be touched upon here, except to 
note that so far the only case with 2-hody interactions in 
which B-E condensation has been rigorously established 
is hard core lattice bosons, but only at half filling [p^ . 

Dyson effectively converted region 1 into region 2. 
We shall make use of his important idea, which substi- 
tutes a very soft potential for the original one (even a 
hard core) at the price of sacrificing the kinetic energy. 

We assume that the N particles are vci a. L x L x L 
cubic box, Vl. The particle density is then p — NL^^. 
It is well known that the energy per particle in the ther- 
modynamic limit, eo(jo), does not depend on the details 
of (reasonable) il, so we are free to use a cube and take 
N —> oo through any sequence we please, as far as eo{p) 
is concerned. We set N — kM with k an integer and 
M the cube of an integer, because we shall want to di- 
vide up f2 into M smaller cubes (called cells) of length 



^ = {k/pf/^. We will take M ^ oo with £ and fc = pt^ 
fixed, but large. 

The N-body Schrodinger operator is 

N 

H=-^lY,^^+ Yl (4) 

i=l l<i<j<N 

For boundary conditions wc impose Neumann (zero 
derivative) boundary conditions on fi. It is well-known 
that Neumann boundary conditions give the lowest pos- 
sible g.s.e. for H, and hence its use is appropriate for 
a discussion of a lower bound for the g.s.e. Denote this 
Neumann g.s.e. hyEo{N,L). 

Now divide into M cells and impose Neumann con- 
ditions on each cell, which, as stated before, lowers the 
energy further. We also neglect the interaction between 
particles in different cells; this, too, can only lower the 
energy because v >0. 

A lower bound for Eo{N, L) is obtained by distributing 
the N particles in the M cells and then finding a lower 
bound for the energy in these cells, which are now in- 
dependent. We then add these M energies. Finally, we 
minimize the total energy over all choices of the particle 
number in each cell (subject to the total number being 
N). Despite the independence of the cells, the latter 
problem is not easy. In particular, something has to be 
invoked to make sure that we do not end up with some 
cells having too large a number of particles and some cells 
having too few. 

With L,N and M = N/{p£^) fixed, let Mc„, for n = 
0, 1, 2, . . . denote the number of cells containing exactly 
n particles. Then the particle number and cell number 
constraints are 

^ nc„ = fc = pf, ^ c„ = 1 (5) 

n>0 n>0 

and our energy bound is 

Eo{N,L) > M min ^ c^ Eo{n,e), (6) 

n>0 

where the minimum is over all c„ > satisfying (|^). 

The minimization would be easy if we knew that 
EQ{n,£) (or a good lower bound for it) is convex in n, 
for then the optimum would be c„ = Sn,k- This convex- 
ity is very plausible, but we cannot prove it (except in 
the thermodynamic limit, where it amounts to thermo- 
dynamic stability). What we do know instead is super- 
additivity: 

Eoin + n'J)>Eo{n,e)+Eo{n',i) (7) 

for all n, n' , and this turns out to be an adequate sub- 
stitute for controlling the large n terms in (|6|). 

Eq. is an immediate consequence of the positiv- 
ity of the potential and it is used as follows. Suppose, 
provisionally, that we have a lower bound of the form 
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Eq (n, t) >K{e)n{n-l) for < n < 4fc, (8) 

with K{i) independent of n for < n < 4fc. In fact, 
we shall later prove that for small enough p (and hence 
small enough k) and suitable £, (||) holds with 

4™r3(i-c" yi/"), (9) 

with C" some constant. [However, the analysis we give 
now, leading to (p^), does not depend on this particular 
form oiK{i).] 

Split the sum in (^) into two pieces: < n < 4fc 
and 4fc < n. Let t = Y^^KAk — that 
k ~ t — X]n>4fc From (g) and Cauchy's inequality 
(and E„<4fc cn < 1) 

CnEo{n,£)>K{£)t{t-l). (10) 

n<4fc 

On the other hand, if n > 4fc then, by (0), Ei^{n,l) > 
{n/8k)Eo{AkJ), so 

J2 ^nE,{n,£) > ^K{l){Ak - 1). (11) 

n>4fc 

Upon adding (|l^ and (|l|) the factor t(t-l) + (fc-t)(2fc- 
1/2) is obtained. Although the number t is unknown, we 
note that this factor is monotone decreasing in t in the 
interval < t < k { which is where t lies, by Thus, 
we can set t — k and obtain the same bound as if we had 
convexity, i.e., 

E^{N,L)>N K{1) {pf (12) 

In summary, if we can show (^ for a box of a fixed size 
for all particle numbers up to n = 4p£'^, then we will 
have obtained our goal, (H), in the thermodynamic limit 
provided we can show that the K in (^ satisfies (^ with 
the constant C" when £ is large compared to the mean 
particle spacing, i.e., p £^ > C" Y^^/^'^ . Then the C in 
(|) is equals C + C" . 

We now focus on a single cell and denote the n coordi- 
nates (xi, x„) collectively by X. The first step in prov- 
ing (|^) is to replace the total potential, X]i<j" ^{"^i ~ ^j)j 
by a lesser quantity, Wv (X) , the nearest neighbor poten- 
tial defined by 

1 " 

i=l 

where j{i) is the nearest neighbor to particle i in the 
configuration X. I.e., particle i 'feels' only its nearest 
neighbor. Hence, we replace H by the smaller operator 

H = T + Wv<H, (14) 

where T = — /i^A; is the kinetic energy in (^. Since 
t) > 0, the g.s.e. of H satisfies Eoirijl) < Eo{n,£). 



To get into the kinetic energy dominated region, we 
wish to replace v in ( p^ ) by a gentler potential U. To 
this end we generalize Lemma 1 of and simplify its 
proof. 

LEMMA 1: Let v{r) >0 and v{r) ^ for r > Rq. Let 
U(r) > 6e any function satisfying JU{r)r^dr < 1 and 
U(r) = for r < Rq. Let B C R'^ be star-shaped (convex 
suffices) with respect to 0. Then, for all functions (j), 

j^^\V4>{^)\'' + [\v(r) - paU{r)m^)\'d^^ > 0. (15) 

Proof: Actually, ( |l5|) holds with /i|V0(x)|'^ replaced 
by the (smaller) radial kinetic energy, u\d(j){x.)/dr\'^ , and 
thus it suffices to prove the analog of (|l^) for the integral 
along each radial line, and to assume that 0(x) = u{r)/r 
along this line, with u(0) = 0. Let us first prove ( |l5| ) 
when J7 is a delta- function at some radius R > Rq, i.e., 
U{r) = R~^S{r — R). Then, it is enough to show, for all 
u, that 

/ p\u'{r) - u(r)/rp + -v{r)\u{r)\^ dr > pa\u{R)\'^ R'^ . 
Jo 2 

(16) 

If the length of the radial line is less than R then ( |l^ ) 
is trivial. Otherwise, normalize u by u{R) = R — a, 
and ask for the minimum of the left side of (^) under 
the condition that u{Q) — 0, u{R) = R — a. This is a 
simple problem in the calculus of variations and leads 
to the scattering length equation (^. If we substitute 
the solution into (p6[), integrate by parts, and note that 
UQ{r) — r—a for r > Rq, we find that ( p6| ) is true \ia < R, 
which is true since uq > Q. Finally, by linearity and the 
fact that U{r) = J r~'^S{r — s)U{s)s'^ds, the (5-function 
case implies the general case. Q.E.D. 

We select our U by picking some R^ Rq and setting 

U{r) = 3(i?3 - Rl)-^ for i?o < r < i? (17) 

and U{r) =0 otherwise. Later on we shall choose R, and 
we shall take 

i?o < < P"^/^ < ^- (18) 

By further decomposing a cube into Voronoi cells 
(which are always convex), Dyson Q| deduces from 
Lemma 1 that H is bounded below by a nearest neighbor 
potential, as in (^3|), i.e., 

H >H> pa W[/(X), (19) 

where Wu is as in (p^), with v replaced by U. For the 
hard core case, Dyson estimates the minimum (over all 
X) of W[/(X), for a U similar to (p7|), and gets a lower 
bound for all p, but 14 times smaller than (|l|). We follow 
another route. An important quantity for us will be the 
average value of Wu (X) in a cell, denoted by ( Wu ) • 
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To compute ( i7(xi— Xj(i)) ), for example it is easiest to 
do the X2, . . . , x„ integrations over the ceU first and then 
the Xi integration. Provided Xi is in the smaher cube 
which is a distance R from the ceh boundary (whose vol- 
ume is {£—2R)^), the probability that Rq < |xj — Xi| < R 
is 47r(i?^ — Rq)/3£^. Thus, performing the Xi integration 
over the smaller cube, and then adding similar contribu- 
tions from C/(x2 — Xj(2))j etc., and using (p^, we get 



{Wu)> 



3n{£ - 2Rf 
(i?3 - i?3)£3 



At: , J 2R 



1 



1 + Q(n-1)' 



with 



Q = 47r(i?3 - Rl)/U^ < 1. 



(20) 
(21) 

(22) 



In (|21|) we used < [l + {n-l)x\-^ forO < x < 1. 

Note that (|l]) is of the form (1). 

By similar reasoning, we obtain the upper bound 



( Wc/ ) < 



3n 



R^-Rl 



in 

Since [/(r)^ = 4:TT{Q£^)^'^U{r), we also obtain 
{WIj)< Ann (Qe^^i Wu ). 



(23) 
(24) 

(25) 



We can now use Lemma 1 and these averages to obtain 
(^ and (H). Instead of using (|l^) alone, we pick some 
< e <C 1 and, borrowing a bit of kinetic energy, define 



H = eT+{l-e) /iaWc/(X) 
By Lemma 1 and u > 0, we have 
H > H > H 



(26) 
(27) 



We shall derive (||) and (^) from a lower bound to H. 

Although e is small, we regard Ho = eT as our un- 
perturbed Hamiltonian and V = {1 — e)/ia>V[/(X) as a 
perturbation of Hq. The ground state wave function for 
Hq is *o(X) = ^-3"/2 and Hq^q = Ao*o = (Neumann 
conditions). The second eigenvalue of Hq is Ai = e[in jP' . 
Note that the ground state expectation, (^'ojWc/l^'o), is 
precisely the average ( Wc/ ) mentioned in (po|)-(p5|). 

Temple's inequality states that when a perturba- 
tion V is non-negative (as here) and when Ai — Aq > 
(^o|V"|5'o) then the g.s.e., i?o, of the perturbed Hamilto- 
nian H = Hq + V satisfies 



En > An 



(*o|^|*o)- 



{^o\V^\^o)-{^o\V\^oy 
Ai - Ao - (*o|l^|*o) 



(28) 



We apply this to our case with Ai — Aq = e^n/£'^ 
and V = {1 — e)fiaWu- We neglect the (positive) term 



(^olVI^-o)^ in (H) and we use (HJ), (||) and (||). We 
also use 1 — e < 1 in two appropriate places and find 



EojnJ) 



>(1-^) 1 



Anan 



1 



en ■ 



(29) 



Apart from some higher order errors, (|29|) is just what 
we need in (||), (||). Let us denote the order of the main 
error by F", and we would like to show that a = 1/17 
sufhces. The errors are the following: 

From the (1 — e) factor, we need e < 0(F"). 

From the Q{n+1) error in (|l|) we need Qp£^ < OiY"). 

From the R/£ error in ^) we need R/£<0{Y°'). 

From (11) we need a£'^{ Wu )/e < 0{Y°') and 
p£^a/{R^e) < 0(r"). 

All these desiderata can be met with e = F", R/£ = 
Y°',Q = 0{Y"), pR^ = anda = 1/17— as claimed. 
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